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Abstract 

We perform fits of Aqcd and the gluon fragmentation function D{x,Q) at initial scale Qq ^ 
Aqcd to charged light hadron momentum spectra data by evolving in the Modified Leading Log- 
arithm Approximation. Without additional assumptions, we achieve a good description of the 
available data for ^ = ln(l/2;) up to and around the Gaussian peak, and values of Aqcd acceptably 
close to those in the literature. In particular, we find that this procedure describes the position of 
the peak, and, in contrast to the Limiting Spectrum, also the normalization. 
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I. INTRODUCTION 



Cross sections in which hadrons are detected in the final state currently cannot be reliably 
calculated from first principles in Quantum Chromodynamics (QCD). However, as a result of 
the factorization theorem, one can separate these cross sections into perturbatively calculable 
hard parts convoluted with parton densities if there are hadrons in the initial state and 
fragmentation functions (FF's), which contain all the information on the soft transition from 
a parton a to the produced hadron h. FF's for charged particles have been well determined 
over large and intermediate values of the hadronic momentum fraction x = 2p/ ^/s, where p 
is the momentum of the hadron h and y/s is the centre-of-mass energy, by fitting to a wealth 
of experimental data However, data at x < 0.1 have always been excluded from fits 
because the convergence of the fixed order perturbation series for the evolution of the FF's 
is spoilt by terms of the form a"ln™'(l/x), and so FF's are not well understood at small 
X. A theory which resums these logarithms at leading and sub-leading order exists — the 
Modified Leading Logarithm Approximation (MLLA) [2] (for reviews see Refs. 0,0). The 
MLLA is a systematic improvement over an earlier approximation, the Double Logarithmic 
Approximation (DLA), which resums leading logarithms by summing tree level diagrams in 
which the outgoing gluons are strongly ordered in their angles of emission, thereby giving 
the largest logarithm of the gluon FF at the order in of the diagram . 

The MLLA has been primarily studied in the context of the Local Parton-Hadron Duality 
(LPHD) approach j^. Here, one assumes that, when the longitudal momentum fraction z 
of the observed hadron relative to the parent parton is low, a sufficiently inclusive hadronic 
process has similar properties to the corresponding process involving partons with transverse 
momentum less than the order of the hadron's mass. The FF's describe all partons with 
transverse momentum less than the factorization scale Q, so for light hadron production the 
shape in x space of the initial FF's with initial factorization scale Qo = O(Aqcd) will be 
similar to the shape of the probability for a parton to emit a parton, i.e. these FF's are delta 
functions in {1 — z), and only the absolute normalization Kh is undetermined. Using this 
assumption, and fixing Qq = Aqcd, where the MLLA resummed evolution is well behaved, 
leads to the so-called Limiting Spectrum 0, which can make predictions for data at small 
X with just two free parameters to be fitted, Kh and Aqcd- Together with the conventional 
choice Q = y/s/2, this approach has been very successful at describing the ^ = ln(l/a:) 
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dependence of small x data, provided some modifications are made to the MLLA evolved 
normalization: in Ref. an additional component not provided by the MLLA was added, 
whereas in Ref. a different Ky^ was fitted for each value of i/i. 

In this paper, we are interested in studying MLLA evolution without using strong as- 
sumptions about the non-perturbative physics such as the LPHD, or modifying the MLLA 
evolution itself. There are a number of important reasons for this. Firstly, it is interesting to 
determine whether the MLLA can describe the \fs dependence of the overall normalization 
of the data. Secondly, in current analyses, where only the NLO calculation has been used, 
such as in Ref. Q], fitting is achievable only to data for which x > 0.1 < 2.3). A continued 
rise in the data as x decreases is predicted, whereas the experimental data reach a peak and 
then fall. Therefore it is important to know if one can use the MLLA to improve the hard 
part at small x such that the fitting can be extended over that in the literature to include 
data for which x < 0.1 and therefore, since the cross section depends on the FF's for 2; > x, 
obtain FF's in that region. Thirdly, using weaker assumptions will allow for a purer test of 
the MLLA and determine its kinematic range of validity better. This can be achieved, as 
in global fits, by taking Qo ^ Aqcd to stay in the perturbative region, in which case one 
does not need to assume the Limiting Spectrum to be valid, and absorbing the soft physics 
at energy scales less than Qo into a parameterized FF, whose free parameters can be fitted 
to data at Q by evolving this initial FF in the MLLA. The distorted Gaussian in £, with no 
MLLA evolution, gives a good description of data over the range of Q ^ Aqcd |Zl , so 
we shall employ this parameterization at Q = Qo- 

The organization of this work is as follows. In Section |n] we shall repeat the basic MLLA 
equation in moment space and discuss various approximations to it. Section IIIII contains 
the comparison with e^e~ single charged hadron spectra at the larger scale Q > Qo and the 
determination of Aqcd- In Section |TVl we make some changes to the theoretical input to 
further understand the limitations of our general approach. Finally in Section |3 we present 
our conclusions. 

II. MLLA EVOLUTION 

Before we present our results concerning the evolution of the low x spectra based on the 
MLLA evolution equations, we shall list the basic equations on which our analysis rests. 
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We work in the LO approximation where the inclusive cross section for e~^e hX as a 
function of x is related to the FF's -D^ (x, Q) for the transitions q —>■ h and q h hj 

1 da'^ _E,elD^{x,Q)_ 



O-tot dx Y.q(^l 

where is the electoweak charge on quark g, citot = X]g 47rega;^/(3s) is the total cross 
section, and is the cross section for the inclusive production of a hadron h. 

As usual we also use the variable ^. At sufficiently small x, i.e. large ^, the contribution 
to the cross section from the non-singlet sector may be neglected in our approach since the 
non-singlet evolution is free from small x logarithms. Writing each quark FF in the form 

Q) = Q) + D^s,,(x, Q), (2) 

where the singlet D^{x,Q) is defined to be the sum over all quark FF's divided by the 
number of quark fiavours Nf and the D^^ ^{x, Q) are the non-singlets, we therefore see that 
each quark FF in Eq. ([Q) may be replaced by the singlet FF. Furthermore, at small x one 
can make the approximation, good within MLLA accuracy, that the singlet FF is related to 
the gluon FF by 

D|(x,g) = ^Dj(x,Q), (3) 

where Cp = {N^ — 1)/{2N^. Using these approximations in Eq. (^, we find that the cross 
section in the MLLA can be written 

(Ttot dx Nc " 

In other words, for describing the fragmentation q{q) h at large ^ we can just use the 
FF for g ^ h. Note therefore that the cross section can only depend on Nf through the 
evolution of the gluon FF, which we will consider just now. In the following we shall skip the 
upper and lower indices and write Dg{x, Q) = D{x, Q). The MLLA equation for D{x, Q) is 
most easily written by introducing the moment transform Dj{Q) of D{x, Q), which is 

DjiQ)= [ dxx^-'Dix,Q), (5) 

^0 



with the inverse transformation 



D{x,Q)= I -L:,-Wj{Q), (6) 



where r must be chosen such that the integration contour hes to the right of all poles in 
Dj{Q). We introduce = j - 1 and write D^{Y) = Dj{Q) with Y = \n{Q/Qo). Then the 
MLLA equation for D^{Y) is Q 

where a = llNc/3+2Nf /{3N^). The solution to this equation for D^{Y) is weakly dependent 
on Nf. Indeed, as shown in Ref. the moments of the data calculated with Nf = 3 
and those calculated with Nj increasing by unity whenever ^/s is large enough for the 
contribution from heavy quark flavour to become relevant give similar results up to ^/s = 
202 GeV within the error range on the moments extracted from the experimental data. 
This observation is also substantiated by a recent experimental analysis ll|, where it was 
found that at the Z° resonance, where the effect of heavy quark production is maximal, 
the ^ spectra at the peak determined for all flavours differs from the one for just the light 
flavours by about 8%. In analyses using the Limiting Spectrum it has been sufficient for all 
available data to set Nf = 3, and we will therefore use this value throughout this paper. By 
introducing the anomalous dimension 7aj(«s), we have 

-Y 



D^{Y) = D^iO)exp(^J^ 



dyiu{as{y)) \. (8) 



If -Dtj(O) is known from the FF at the starting scale Qq, which must be taken from experi- 
mental data, Eq. (jSJ gives us the solution for arbitrary F, if we know 7^(as). Equation ((Tj) 
is equivalent to the following differential equation for 7^: 

{uJ + 7..) lu. - 4iV,^ = _ a + 7^) ^ + a6 (^)\ (9) 

ITT das 27r \27r/ 

where 

d 

= = -6^, (10) 

with h = llA''c/3 — 2Nf/3. The first term on the right hand side of Eq. originates from 
the running of a^. The second term gives the hard single-logarithmic correction to the DLA 
soft emission. The last term is formally a next-to-MLLA term, which may be neglected. A 
general solution of Eq. ((Tj) in terms of confluent hypergeometric functions is known js, Q|. 
Equivalently, one can solve Eq. Q in terms of Whittaker functions. However, since the 
MLLA equation is only valid in the region <^ 1 and u = 0(^/aJ), we can obtain a 
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simpler but equally accurate solution to Eq. Q by expanding in asjoj <^ 1 while keeping 
as/uj^ = 0{l) fixed, 



n=l 



and solving for each term. The first and second term will then resum double and single 
logarithms respectively, but the higher terms obtained this way will be incomplete since the 
MLLA does not treat terms in which are of 0{q^J'^) or higher in the region of validity 
given above. 

The DLA corresponds to the n = 1 term only in Eq. in which case the terms in Eq. 
© proportional to (3{as) and a can be neglected. One obtains two solutions: 

1 



7^ = 2 ( ± + H ] , (12) 



with 



ll = 4iV.g. (13) 
7^ = = 7^ = -2w, (14) 

UJ u Ztt 

i.e. 7^^ has the familiar singularity ~ l/cj which determines the small x behaviour of D{x, Q) 
in the Leading Logarithm Approximation (LLA). Therefore the correct solution in the DLA 
is 7(^ = 7^- This solution is finite for — ^ and is equal to 70 ~ ^/o^s■ 

Once the solution for the n = 1 term in Eq. |TT| has been chosen, there is only one 
solution for the n = 2 term, and we have finally 



For a, — > we obtain 



7. = i(-- + ^^)+| 



7o « / 1 , ^ 



1 + 



UJ J Uj"^ 



+ 472 2 \^ ^^pTa^I 

This approximate solution is usually referred to as the MLLA result p|. The term propor- 
tional to a modifies the — > limit to 



which reproduces the finite correction to the LO 7+ in the LLA. The result in Eq. (jl5j) must 
be substituted in Eq. (jH)) to obtain the corresponding MLLA solution for Di^{Y). Writing 

D^{Y) = D^{0)D^{Y), (17) 
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we have 

lnD^{Y)= [ dyj^My)). (18) 

Using the LO formula as{y) = 2TT/[b{y + A)], where we introduce A = ^(Qo/Aqcd), the 
integration in Eq. (fTSj) with given in Eq. (fT3j) yields 

lnDUY) = fiu,Y,X)-fiu,0,\), (19) 

where 

fiu, Y,\)=-^Z+ ^^/Z{Z + 4A) + {2A - B) In (VZ + AA + Tz) 

fl B\ 1 (20) 

+ (4" 2") ^^Z--HZ + AA). 

In Eq. fOUj) we introduced A = ANc/ihoj), B = a/b and Z = a;(F + A). Then the solution 

Duj(Y) can be written as 

D^{Y) = e/(-'^'^)-/(-'0'^), (21) 

with 



VZ + AA + ^/Z 



.1 2A-B 



s 



Z-^. (22) 



By fixing Aqcd the evolution of D^(Y) is completely determined by Eqs. (PT|) and 
This solution has for Y ^ 00 the following asymptotic behaviour: 

g/(.,y,A) ^ ^A-B_ (23) 

In Refs. S and Q it was found that for > 1, 7^ in the MLLA accidentally mimics 
the behaviour of the LLA LO 7^^ reasonably well. This is aided by the observation that the 
cij — > 00 limit of Eq. (fT3j) is equal, up to terms of 0{1/lj), to that of Eq. (fTBj) . the as, — > 
limit of the LLA LO 7^^, whose 0{uj) corrections turn out to be rather unimportant at 
uj = 0(1) and negative at large u. Therefore we neglect those corrections beyond MLLA 
which are important at small ^. 

Solving Eq. Q for the n = 3 term of Eq. ()11|) gives us a contribution to the next-to-MLLA 
correction which reads 

ao\3 /a. 



UJ J \UJ 

a 



2tx 



2 7o ah I 1 uj^ 



The addition of this term to the expression in Eq. ()15|) would give a more accurate ap- 
proximation to the exact solution to Eq. Q if Eq. (fTT|) were a suitably convergent series. 
However, these results are not complete at next-to-MLLA order; in particular they refer only 
to gluon jets. In the complete next-to-MLLA cross section both the evolution and Eq. 
obtain a correction, the latter arising from the energy dependent differences between quark 



and gluon jets 13|. In any case, we can at least use the term in Eq. ()24|) to determine the 
stability of our form for the ML LA evolution in different regions of C, and ^/s. 

Finally, since partons are treated as massless in the MLLA, the parton momentum spec- 
trum is equivalent to the parton energy spectrum. Consequently the MLLA formalism needs 
modification in order to incorporate hadron mass effects, which become more relevant as 
^ increases. However, such effects will be neglected in our analysis since otherwise model 
assumptions are needed. 

III. FITTING THE EXPERIMENTAL DATA 

In this section we test how well the MLLA described in the previous section agrees with 
experimental data, both by fitting free parameters to data sets and by using the resulting 
fitted parameters to predict other data sets. 

Aqcd is the only parameter on which MLLA evolution depends, and should therefore be 
obtainable by fit ting to data at widely separated energies, starting from TASSO data at 
Qq = 14 GeV/2 IJ]. Therefore we use data at the highest namely the recent data at 
^/s = 202 GeV from OPAL [9], as well as data at 91 GeV 15| from the same collaboration, 
which have the highest accuracy. For all experimental data used in this paper, systematic 
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we impose a lower bound on 



and statistical errors are added in quadrature. As in Ref. 
the OPAL data of 

e > 0.75 + 0.33 In (v^/GeV) (25) 

since at lower ^ the experimental errors are too small to fit using only one parameter. 
Including these small ^ points in fact does not change the results significantly but leads to 
a much higher minimized x^- To control the number of data used in the non-perturbative 
region of hadronic momentum p = O(Aqcd); we introduce a cut-off mass scale m and impose 
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an upper limit on the data used of p > m, or 

f<l.>^. (26) 

The initial gluon FF used for this fit was obtained by independently fitting it to data at the 
lowest ^/s, namely the TASSO data at y/s = 14 GeV, using a distorted Gaussian, 



xD{x,Qo) = — ^=exp 
a\/2n 



8 2 4^ ^ 6 24 



(27) 



where 5 = [C, — ^)/cr and Qo = 14 GeV/2, and the results are shown in Table HI The errors 
were obtained from the diagonal components of the inverted matrix of second derivatives of 
at the minimum. Since this method assumes that is quadratic in the parameters, these 
errors should not be taken too seriously. In this case there was no need to impose a lower ^ 
bound on the TASSO data since there were 5 free parameters in the fit. We also imposed no 
upper C, bound on this data, since doing so either made little difference for m < 0.5 GeV or 
did not constrain the parameters sufficiently for m > 0.5 GeV. The achieved per degree of 
freedom, Xdf' is 0.76, and the results in Table IJ for the parameters of the distorted Gaussian 
fit agree well with earlier fits in the literature |l6l |. 

TABLE I: Fit of a distorted Gaussian to all 20 TASSO data points at ^/s = U GeV with Qq = 



Parameter N ^ a'^ s k 

Value 9.71 2.33 0.61 -0.11 -0.77 
Error 0.11 0.01 0.02 0.05 0.12 



The resulting values of Aqcd when performing this procedure, and cutting the data using 
values of m ranging from 0.3 to 0.6 GeV, are shown in Table ITU where it can be seen that 
the obtained value for Aqcd with Nf = 3 depends somewhat on the upper limit for C,. The 
errors were calculated by varying Aqcd, in both directions, from its value at the minimum 
until increased by unity. The errors were found to be symmetric and close to the inverse 
of the second derivative of with respect to Aqcd- If we choose the Aqcd from the fit with 
the smallest Xdf have Aqcd = 317 MeV in reasonable agreement with LO Aqcd values 
with Nf = 3 obtained in other analyses 
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TABLE II: Four independent fits of Aqcd to OPAL data at 91 and 202 GeV, where the cuts in 
each case are labelled by the value of m. 



m (GeV) 0.3 


0.4 


0.5 


0.6 


Aqcd (MeV) 258±8 


293±9 


307±10 


317±10 


xIf 7.0 


3.0 


2.3 


1.8 



The fits for m = 0.3 and 0.6 GeV are shown graphically in Figs, ^andj^ These figures 
also show the predictions of the respective fits for TPC data at 29 GeV [l2|, TASSO data 
at 35 and 44 GeV Q, TOPAZ data at 58 GeV and OPAL data at 133 3 and 172 
GeV [20]. In all plots in this paper, each curve is shifted up from the curve below by 0.8 for 
clarity. The data are well described almost up to the peak, about one half or one unit in ^ 
below. Beyond the peak the predictions fail. 

In Ref. |8| a much better agreement with data over the whole ^ range was obtained, but the 
MLLA prediction was modified in two aspects. Firstly, an energy independent background 
term was added to the MLLA multiplicity formula and, secondly, a correction for mass effects 
at large ^ was added. In contrast, our fit, starting from the Gaussian parameterization of the 
TASSO data at 14 GeV, predicts the ^ distributions at higher energies using only a single 
parameter, Aqcd! it is remarkable that MLLA evolution predicts all data sets at higher 
energies very well up to the peak region. The discrepancy beyond the peak region in our 
case is too large to be attributed to mass effects. 

The distorted Gaussian parameters are highly correlated with one another, so in order to 
constrain them better, we fit these parameters and Aqcd simultaneously to TASSO data at 
14 GeV and OPAL data at 91 and 202 GeV. We again impose the upper bound of Eq. ()26p. 
but this time on all three data sets for consistency. However, we impose no lower bound on 
any of the data, because all distorted Gaussian parameters are free in the fit. Taking first 
m = 0.4 GeV, we obtained the results presented in Table ITTTl with Xdf = 2.3, which are 
shown graphically in Fig. 01 This figure also contains predictions for other data sets not used 
in the fit. The case for m = 0.5 GeV is shown in Table HVl and Fig. 01 where Xdf = 2.1. 
Since the dependence of the parameters cannot be adequately approximated by a 

quadratic, and since the present study does not aim at a precise determination of Aqcd, we 
refrain from calculating the errors in these and subsequent tables. 
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FIG. 1: Fit of Aqcd to OPAL data at ^ = 91 and 202 GeV, after fitting the initial gluon FF 
to TASSO data at 14 GeV. The ^ region of data is chosen as described in Eqs. H25() and 
and is indicated by the vertical dotted lines. The upper bound corresponds to m = 0.3 GeV. The 
predictions from this fit of other data sets is also shown. The lowest curve shows the independent 
distorted Gaussian fit to TASSO data at 14 GeV. Each curve is shifted up by 0.8 for clarity. 

TABLE III: Fit of gluon FF and Aqcd to all TASSO data at 14 GeV and OPAL data at 91 and 
202 GeV (88 data points), with m = 0.4 GeV. 

N I s k Aqcd (MeV) 

7.86 2.11 0.40 -0.46 -1.32 649 



TABLE IV: As in Table |TTT| but with 83 data points and m = 0.5 GeV. 

N I s k Aqcd (MeV) 

11.80 2.60 0.67 -0.26 -1.48 87 



The data around the peak region are better described for m = 0.5 GeV. We note that 
there is a large difference between the parameters in each case, which may be due to the fact 
that the theory cannot accomodate some or all of the three main features of the data, being 
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FIG. 2: As in Fig. ^ with an upper bound in on the data used corresponding to m = 0.6 GeV 
and indicated by a vertical dotted hne. Each curve is shifted up by 0.8 for clarity. (Note again 
that the lowest curve is from an independent fit to all TASSO data at ^/s = 14 GeV, and hence is 
identical to the corresponding curve in Fig. 



the position of the maximum, the width and the normahzation. Indeed, two local minima 
were found in each fit, and the global minimum shifted from one of these two local minima 
to the other as m increased from 0.4 to 0.5 GeV. 

The resulting values for Aqcd in the fits of Tables lllll and IIVI are clearly too different 
for either of them to be taken seriously. This is probably due to the fact that Aqcd and 
the distorted Gaussian parameters are highly correlated with one another. In the fits of 
Table ITTl the values of Aqcd are more consistent with each other since they were completely 
uncorrelated with the distorted Gaussian parameters. 

To better constrain all the parameters would require using more data sets in the fit. 



Therefore we fit the distorted Gaussian parameters and Aqcd to all avai 



namely the data sets in Figs. 



DELPHI M, L3 



as well as TASSO data at 22 GeV |2l|, ALEPH 



able data sets, 



2JI and SLD "25] data at 91 GeV, ALEPH data at 133 GeV [261 , DELPHI 
data at 161 GeV [271 and OPAL data at 183 and 189 GeV Q- For m = 0.5 GeV, we 
obtained the results shown in Table IVl and Fig. for which Xdf — ^-0 was achieved. The 



12 




1 2 3 4 5 6 

FIG. 3: Fit of gluon FF and Aqcd to TASSO data at 14 GeV and OPAL data at 91 and 202 
GeV, with an upper bound in ^ on the data used corresponding to m = 0.4 GeV and indicated by 
a vertical dotted hne. Other data sets are shown for comparison. The upper bound on ^ for each 
data set used in the fit is indicated by a vertical dotted line. 

results do not differ significantly from those in Table IIVI and Fig. El nor from similar fits 
with m = 0.4 and 0.6 GeV, for which we obtained Aqcd = 106 and 129 MeV respectively. 
In all cases we found that there were more than one local minimum, from which we selected 
the minimum with the smallest x^- 

TABLE V: Fit to all available data (413 data points), with m = 0.5 GeV (see text). 
N I s k Aqcd (MeV) 

11.65 2.57 0.70 -0.19 -1.17 130 



The fit in Fig. El is the main result of this paper. A global fit in which the parameters of 
the distribution at the lowest scale Qq are fitted simultaneously with the parameter Aqcd 
leads to an improvement over the fits in Figs. Q and El At all energies the description is 
now good up to the peak or even beyond. We stress again that this fit, beyond the input 
parameterization at Qq, does not involve any additional assumptions beyond the MLLA 
evolution. 
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1 2 3 4 5 6 

FIG. 4: As in Fig. El but with m = 0.5 GeV. 
IV. FURTHER STUDIES 

In all our fits so far we obtained a good description of the data below the maximum, i.e. 
for small ^ values, but a rather bad description in the region above the peak, i.e. for large ^ 
(with the exception of the TASSO 14 GeV data, when it was fitted over the whole ^ range). 
This discrepancy may have several reasons. 

Since the MLLA approach is supposed to be particularly valid for sufficiently large ^, 
presumably in the peak region, despite the discussion at the end of Section [TTl it may have 
been necessary to exclude data below a given ^, e.g. that of Eq. in our approach of 

fitting the distorted Gaussian parameters and Aqcd simultaneously to all three data sets. 
However, with this approach only a few data points are left, in particular for the TASSO 
data at 14 GeV, when imposing also the upper limit on ^ in Eq. fl26j) . Therefore a lower ^ 
cut with the approach applied in Tables ITTTl and HVI does not work. 

Alternatively, it may be that the upward evolution of the higher moments tends to become 
unstable. To investigate this possibility, we fit the distorted Gaussian parameters and Aqcd 
to TASSO data at 14 GeV and OPAL data at 91 and 202 GeV as before, but this time 
we set Qq = 202 GeV/2, i.e. we fit the initial distribution at the highest energy and evolve 
downwards. The results of this fit are shown in Table IVll and Fig. El where Xdf — 3-9. The 
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Q I I \ I \ I \ I \ I - I I 

1 2 3 4 5 6 

FIG. 5: Global fit of gluon FF and AqcD) with m = 0.5 GeV, to the data shown here and other 
data listed in the text. 

value of Aqcd obtained is in good agreement with that obtained in other analyses. The 
resulting ^ distributions also fit better to the data at the highest energies at larger ^ values 
beyond the peak, whereas at the remaining energies the description of the data around the 
peak becomes worse. 

TABLE VI: Fit of gluon FF to TASSO data at 14 GeV and OPAL data at 91 and 202 GeV (83 
data points), using Qq = 202 GeV/2 and downward evolution, with m = 0.5 GeV. 

N I s k Aqcd (MeV) 

26.83 3.66 1.17 -0.52 -1.49 225 



Another possibility for our large discrepancy may be due to the region of function space 
in ^ available to the parameterization in Eq. ()27j) being insufficient. To enlarge this region, 
we added a term 

C5(5' + Cq5^ (28) 

to the argument of the exponential in Eq. ()27|) , and include C5 and Cg in the list of parameters 
to be fitted. However, when performing this fit with m = 0.5 GeV, there was no significant 
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L^.ti 1 , 1 , 1 , 1 ^ ^J^Ju — ^ 

1 2 3 4 5 6 

FIG. 6: Fit of gluon FF to TASSO data at 14 GeV and OPAL data at 91 and 202 GeV, using 
Qo = 202 GeV/2 and downward evolution, with m = 0.5 GeV. 

improvement over the fit in Fig. H 

This indicates that the failure to describe the region above the peak is inherent to the 
MLLA formahsm as it is apphed here. A better approximation to the full analytic solution 
to Eq. ((7j) may improve the large ^ description, since it includes certain corrections of next- 
to-MLLA order. Such corrections are also included, for example, in the Limiting Spectrum 
within the LPHD approach, where, compared to our fits, a better description of the data 
beyond the peak is achieved, given suitable modifications to the MLLA evolution of the 
normalization (see Section Therefore we repeated the fit of Fig. ^ but this time including 
the extra term given by Eq. (j24|) in the evolution. In this case Xdf increased to 2.6, and this 
increase can be attributed to the fact that the deviations from the data were slightly larger 
beyond the peak. However, up to the peak the description was as good as the fit of Fig. |3 
Furthermore the theoretical curves were rather similar to those of Fig. 0] in the ^ range of 
the data. This suggests that the MLLA can only describe data up to the peak, and that 
a full next-to-MLLA calculation is required beyond the peak, which includes, in particular, 
the correlation between the evolution of quark and gluon jets. 
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V. CONCLUSIONS 



In this work we perform fits to the available momentum spectra data of e^e~ annihilation 
in the energy range 14 - 202 GeV using MLLA evolution. No additional assumptions, such as 
the LPHD, is used other than a conjectured functional form for the gluon FF, and therefore 
we have achieved a particularly pure test of the MLLA. We find a good description of the 
data in the region up to the maximum of the distribution in the scaling variable ^, with 
only a minimal number of parameters. In particular we find that MLLA evolution without 
additional input gives a good description of the normalization up to the peak, and also the 
approximate position of the peak. 

Our fitted values of Aqcd cover a large range. However, in our model-independent ap- 
proach, there is some theoretical ambiguity in Aqcd- We have chosen the renormalisation 
and factorization scales to be Q = y/s/2, but we could also have chosen some factor of this, 
of 0(1). With this theoretical error, our results for Aqcd are consistent with those of other 
studies 

Clearly, our form for the MLLA evolution is insufficient to describe the data above the 
peak. The inclusion of the next-to-MLLA contribution, Eq. (j23)) did not improve our results. 
At this order, a full treatment of momentum distributions would include quark-gluon mixing, 
which may be the most important effect at this order and therefore may significantly help 
to reduce the large ^ discrepancy. 

Finally, it will be interesting to incorporate the MLLA into the full NLO fits which apply 
to the large x range, in order to extend the region of validity towards lower values of x. Our 
recipe for fitting the fragmentation functions is consistent and compatible with the standard 
fitting. 
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